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CONSTRAINTS ON HELIUM ENHANCEMENT IN THE GLOBULAR CLUSTER M3 (NGC 5272): 

THE HORIZONTAL BRANCH TEST 
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ABSTRACT 

It has recently been suggested that the presence of muhiple populations showing various amounts of helium 
enhancement is the rule, rather than the exception, among globular star clusters. An important prediction of 
this helium enhancement scenario is that the helium-enhanced blue horizontal branch (HB) stars should be 
brighter than the red HB stars which are not helium-enhanced. In this Letter, we test this prediction in the 
case of the Galactic globular cluster M3 (NGC 5272), for which the helium-enhancement scenario predicts 
helium enhancements of > 0.02 in virtually all blue HB stars. Using high-precision Stromgren photometry and 
spectroscopic gravities for blue HB stars, we find that any helium enhancement among most of the cluster's 
blue HB stars is very likely less than 0.01, thus ruling out the much higher helium enhancements that have been 
proposed in the literature. 

Subject headings: stars: abundances — Hertzsprung-Russell diagram — stars: evolution — stars: horizontal- 
branch — (Galaxy:) globular clusters: general — (Galaxy:) globular clusters: individual 
(M3 = NGC 5272, M13 = NGC 6205) 



1. INTRODUCTION 

Globular star clusters (GC's) have traditionally been as- 
sumed to be excellent approximations to so-called "simple 
stellar populations," which are idealized systems in which all 
stars were formed at precisely the same time, from a chem- 
ically homogeneous cloud. However, recent observations, 
both photometric and spectroscopic, have cast serious doubts 
on this long-standing paradigm. 

The presence of large abundance anomalies was first iden- 
tified in the GC w Centauri (NGC 5139). In this clus- 
ter, not only such light elements as C, N, O, F, Na, Mg, 
and Al, but also the Fe-peak, s-process, and r-process ele- 
ments, are se en to vary by large amounts from one star to 
the next (e.g., iJohnson et alj |2008, and the extensive list of 
references provided therein). Many of these abundance pat- 
terns seem to extend all th e way down to the main sequence 
(e.g.. [Stanford et aT]|2008l) . This strongly suggests that such 
abundance anomalies owe their origin to multiple star for- 
mation episodes within the cluster, each accompanied by a 
corresponding enrichment of the intracluster medium by the 
ejecta of the massive and intermediate-mass stars which were 
formed in the prior stellar generations. While a large spread in 
Fe-peak abundances has only been detected in lo Cen, many of 
uj Cen's abundance anomalies have also been found in other 
clusters, alb eit often at a (much ) less dramatic level (e.g., 
ISneden et al.ii2004; .Johnson et alj|2005; Yong et aL.2009; see 
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also lGratton et alJl200 4). for a recent review). 

Analyses of the multiple main sequences found in the deep 
color-magnitude diagrams (CMD's) of uj Cen strongly sug- 
gest that the helium abundance Y in the cluster may also 
have changed dramatically from one star formation episode to 
the ne xt (e.g.. iNorris„2004j:lD'Antona et al.ll2005l: iPiotto et alj 
120051) . Recent evidence suggests that other clusters may also 
show sizeable Y variations. According to the deep CMD 
analysis of NGC 2808 by Piotto et al. (2007), multiple pop- 
ulations with different helium abundances are also present in 
this cluster Strong arguments have also been raised in fa- 
vor of helium enhancements among at le ast some of the stars 
in NGC 63 88 and NGC 6441 (see, e.g.. iCatelan eTaP 120061: 
ICaloi & D^A ntona 2007, and references therein). Interest- 
ingly, these four GC's figure among the most massive of all 
Galactic GC's. 

Very recently, it has been suggested that multiple star for- 
mation episodes in GC's, accompanied by widely different 
amounts of helium enric hment, are in fact not th e excep- 
tion, but instead the rule (D' Anto na & Caloill2008h . In the 
particular case of M3 (NGC 5272), such a c l aim ha d previ- 
ously been made also by Caloi & D'Antona (2008), on the 
basis of an analysis of the period distribution of M3's RR 
Lyrae variables and the color extension of the HB blueward 
of the instability strip. Due to the recognized need to as- 
sume a very sharply peaked mass distribution to account for 
the sharply peaked shape of the RR Lyrae period distribution 
inM3 ( Rood & C rocker 1989; Catelan 2004; Castellani et al] 
120051) . Id 'Antona & Caloii suggested that the mass distribu- 
tion of HB stars is always sharply peaked, and that the color 
spread routinely observed in GC CMD's is instead due to in- 
ternal variations in the helium abundance Y. It is important 
to emphasize that such a spread in Y is not required to ex- 
plain the observed period distributions. Rather, the case for 
helium enhancement, at least in the case of M3-like clusters, 
rests almost entirely on the presence of color spreads among 
HB stars. 

If confirmed, this scenario would not only have major im- 
plications for our understanding of how GC's form, but would 



2 



M. Catelan et al. 






Fig. 1 . — Panel a: compaiison between fiducial HB loci (ZAHB, MAHB, 90AHB, TAHB; see definitions in Q for three different metallicities (black lines: 
Z = 0.002; gray lines: Z = 0.001; blue lines: Z = 0.0005), in the My, (fo-y)o diagram. To produce this plot, the Z = 0.001 loci were shifted in My by +0.075 mag, 
and the Z = 0.0005 loci by +0. 140 mag. Note that with these shifts all three of the ZAHB's coincide, except at the extreme red end of the HB, where no unevolved 
stars are found in M3. Panel b: as in panel a, but for a fixed metallicity and two different helium abundances {black lines: Y = 0.23; gray lines: Y = 0.28), and 
without applying any shifts to the CMD positions. Panels c and d: as in a and b, but in the logg, logTeff plane, and without applying any shifts to the data. 



also represent a major shift from the canonical pa radigm, 
dominant since the late-1960's/earlv-197 0's (Castella ni et alJ 
[T969t llben & RoodlfTOTOl: iFaulkneiillQTl iRood, 1 973.) . which 
ascribes the color spread seen among HB stars to the stochas- 
tic nature of mass loss on the red giant branch (RGB). In the 
canonical scenario, blue HB stars lose significantly more mass 
as they climb up the RGB than do their red HB counterparts. 
In contrast, in the .D'Antona & Caloi (.2008) scenario blue HB 
stars, even at the "horizontal" level of the HB, owe their blue 



colors to a higher initial Y. In the specific case of M3, the 
predicted He enhancement, even in the "horizontal" blue HB 
stars, should fa ll in the range between 0.0 2 and 0.025, accord- 
ing to Fig. 4 of lCaloi & D' Antonal JlOOSh . 

Fortunately, the Y enhancement scenario can be tested, us- 
ing a variety of photometric and spectroscopic tools, espe- 
cially for these moderately cool blue HB stars at the "hori- 
zontal" level. In particular, it is well known that HB stars 
with higher Y are brighter, at a given Teff, due to their stronger 
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H- burning shells (e.g.. 



the lD'Antona & Caloi 



Sweigart & Grossll976l) . Therefore, in 



(b-y)o 

Fig. 2. — Comparison between fiducial HB sequences (see text) andtlie em- 
pirical data for M3, in the My, {b-y)Q plane (panel a), Mj,, (h—y)o plane (panel 
b), and My, (fe— v)o plane (panel c). The empirical data were shifted vertically 
so as to produce a satis factor y match to the theoretical red ZAHB, as given by 
the VandenBerg et alj {2003) models for [Fe/H] = -1.61, [a /Fe] = +0.3, and 
Y = 0.236, in the My, {b-y)o plane. The I VandenBerg et al.l RGB sequence is 
also displayed. 
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scenario bluer HB stars should 
be brighter than their redder counterparts. While this may 
help account for the sloping natures of the HB's in the mod- 
erately metal-rich GC's NGC 6388 and NGC 6441 (e.g., 
iBusso et"al]|2007l and references therein), such a CMD test 
has n ever being carried o ut for most of the GC's that, accord- 
ing tol D'Antona & Caloil present multiple He-enhanced pop- 
ulations. 

The purpose of this Letter is accordingly to present a first 
test of the helium enhancement scenario, in the specific case 
of M3, by comparing high-precision photometry for the clus- 
ter in the ,Stromgren (1963) system, as well as spectroscopi- 
cally derived gravities, with theoretical HB models computed 
for a variety of Y values. 

2. OBSERVATIONAL DATA 

The CMD data used in this paper were taken from 
iGrundahl etall (119981 119991) . to which the reader is referred 
for additional details regarding the calibration of our pho- 
tometry. Briefly, our M3 photometry is based on a series of 
images obtained on the Nordic Optical Telescope, using the 
thinned AR-coated 2048 x 2048 pixel HiRAC CCD camera, 
with 0.11" pixel size, thus covering approximately 3.75' on 
a side. Two overlapping fields were observed, with one field 
centered on the cluster center to ensure a large sample of HB 
and RGB stars. 

In addition to the CM D data , we also used the gravities and 
temperatures derived by lBehd (g003), on the basis of observa- 
tions using the HIRES cross-dispersed echelle spectrograph 
on the Keck I telescope. 

3. THEORETICAL MODELS 

In the present paper , we u se th e evolutionary tracks com- 
puted bv lCatelan et all (Il998h and lSweigart & CatelanI (Il998h 
for scaled-solar, heavy-element abundances Z of 0.0005, 
0.001, and 0.002. Evolutionary tracks for a He abundance 
Y = 0.23 by mass were computed for each of these Z values. 
Additional He-enhanced tracks for Y = 0.28, 0.33 were also 
computed for Z = 0.002.^ The th eoretical models were trans- 
formed to the IStromgrenI (11963b uvby system by using the 
color tra nsformation s and bolometric correction tables pro- 
vided by iQemetaD (2004). T he same procedures wer e also 
adopted in the recent work by ICatelan & Co rtds (2008) and 
ICortes & Catelan (2008), where the period-color and period- 
luminosity relations of RR Lyrae stars in the IStromgrenl sys- 
tem were present ed. In addition, w e also u se a set of mod- 
els computed by IVandenBerg et al.l ( 120061) . for a chemical 
composition [Fe/H] = -1.614, [a/Fe] = +0.3 (which trans- 
lates, according to those authors, to a Z = 8.45 x 10""^), and 
7 = 0.236. 

M3 has a metalHcity of [Fe/H] = -1.57 in thelZinn & Wesd 
"1984) scale, and [Fe/H] = -1.34 in the Carretta & Gratto^ 
1997i) scale. Taking into account an enhancement of the a- 
capture elements by [a/Fe] = +0.27 (Carney 1996), and using 
the relation between Z, [Fe/H], and [a/Fe] from Salari s et"al] 
([T99l . we find an overall metallicity ofZ=8.3 x 10 ^ and 
Z = 1.4 X 10"^, respectively, on these two metallicity scales. 
Therefore, our models comfortably bracket the range of pos- 
sible metallicity values for the cluster. 

^ Note that these helium abundances refer to the initial main sequence val- 
ues, and thus do not include the small increase that occurs during the dredge- 
up phase on the RGB. 




Fig. 3. — As in Figure|2] but plotting only the ZAHB's, for the several different Y values as indicated. The thin solid lines indicate interpolated ZAHB loci for 
Y values between 0.23 and 0.28, in intervals of 0.01. Panels d through /present expanded views of the rectangular box around the blue HB "knee" in panels a 
through c, respectively. 
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Since we are primarily interested in constraining the change 
in Y between the red and blue HB, the exact choice of Z 
value is basically irrelevant for our purposes. This is shown in 
Figure [T}?, where our models for Z = 0.0005, Z = 0.001, and 
Z = 0.002 are compared in the M,,, {b-y)o diagram. To pro- 
duce this plot, we first registered the low-metallicity zero-age 
HB (ZAHB) sequences to the Z = 0.002 ZAHB by shifting the 
Z = 0.0005 ZAHB by +0.140 mag and the Z = 0.001 ZAHB 
by +0.075 mag in My. We then applied the same shifts to 
the evolved HB sequences for each Z value, and derived sev- 
eral fiducial loci representing evolved HB stars, as follows: 

i) MAHB, standing for the middle-age HB, or HB ridgeline, 
which gives the average position occupied by all HB stars, as- 
suming a uniform mass distribution along the entire ZAHB; 

ii) 90AHB, or 90%-age HB, which is approximately the lo- 
cus below which one should find « 90% of all HB stars;** 

iii) TAHB, or terminal-age HB, which is simply the He ex- 
haustion locus. 

This registration procedure clearly leads to an excellent 
match over the entire range of ZAHB colors and HB lumi- 
nosities, except (as expected) at the extreme red end of the 
HB and close to the TAHB. Since M3 lacks unevolved HB 
stars at colors redder than about ib-y)o ~ 0.34, and since the 
TAHB lies very far away from where most HB stars are found 
(see Fig. |2]i, these differences are of no concern for our pur- 
poses. Similar results follow when using Mb and M,,, but not 
M,i, which we have found to be affected by strong metallicity 
effects. Since the v- and (especial ly) the u- band bolometric 
corrections from Iciem et al.l (12004) may be quite uncertain 
(D. A. VandenBerg, priv. comm.), we have decided n ot to in- 
clude o ur M-band photometry in this Letter (but see ICatelanI 
l2009bl) . We will, however, include our v-band data, but cau- 
tion that the comparison of these data with the models is less 
reliable than for the yorb data. 

While the metallicity effects are clearly mild, both the 
ZAHB and the evolved HB loci depend strongly on Y, as can 
be easily inferred from Figure [117. This shows that a differen- 
tial comparison between the luminosities of red and blue HB 
stars, around the metallicity of M3, should provide a strong 
indicator of any possible helium enhancement, with only a 
very mild dependence on metallicity. 

Panels c and d in Figure [T] reveal the effects of metallicity 
and Y, respectively, on the logg-logTeff plane. This again 
shows that metallicity plays but a minor role in defining the 
position of a star on this plane, compared with Y. 

4. RESULTS 

We compare the model predictions for fixed (low) Y with 
the empirical CMD data in Figure |2] To produce these plots, 
we have first corrected the empiri cal data for r eddening us- 
ing a EiB-V) = 0.01, taken from Harris] (119961) . and the ex- 
ti nction coefficients for the [Stromgren system summarized 
in ICatelan & CortesI (2008). We then shifted the data verti- 
cally, as required to match the th eoretical red ZAHB com- 
puted by IVandenBerg et al.l ( 120061) to the lower envelope of 
the o bserved distribution in the My, (b-y)o plane. Note that 
these iVandenBerg et al.l evolutionary tracks also provide an 
excellent fit to the RGB of the cluster. This results in a dis- 
tance modulus of (m-M)o = 14.984 for M3. We then applied 

^ In like vein, one may define the "50%-age HB" locus, or 50AHB, as the 
locus occupied by HB stars of different masses which have completed 50% of 
their HB evolution. Note, however, that this is not exactly coincident with the 
MAHB ; since the distribution of HB luminosities is not Gaussian, the MAHB 
is in fact slightly more luminous than the corresponding 50AHB. 
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Fig. 4. — Comparison between predicted and observed loci in the logg, 
logTeff plane, for a metalhcity Z = 0.002 and three different helium abun- 
dance values: 0.23 (panel a), 0.28 (panel b), and 0.33 (panel c). Empirical 
data are provided for M3 (black circles) and Ml 3 (gray circles), from Belj 
l2003j). The blue line in panel a shows the ZAHB locus of iVandenBerg et al] 
12003) models for a F = 0.236, [Fe/H] = -1.61, [a/Fe] = +0.3. 

a shift of Amag = -0. 16 to ourZ = 0.002 models, in order that 
they would match the lVandenBerg et"an ones. As can be seen, 
such a procedure leads to an excellent morphological match 
between the lVandenBerg et al.l ZAHB forZ = 8.45 x lO"'^ and 
our own for Z = 0.002, especially in the My, {b-y)o and M^, 
ib-y)Q planes, except again at the very end of the red HB dis- 
tribution, where no unevolved HB stars are found in M3. 

Except perhaps for a small deviation of the lower envelope 
of the blue HB stars in the immediate vicinity of the "knee" 
from the theoretical (single-K) ZAHB, this figure shows a re- 
markable overall agreement between the model predictions 
and the observations, without an immediately obvious need 
to invoke an increase in Y for the blue HB stars. In addi- 
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tion, all the other HB fiducials for a fixed Y, including the 
MAHB, 90AHB, and TAHB, appear to describe the empir- 
ical data fairly well, without significant evidence for an ex- 
cessive number of overluminous blue HB stars, as would be 
expected in the helium-enhancement scenario (see Fig. (TJ?). 
To be sure, the evolutionary lifetimes along the canonical 
tracks do not provide a per fect match to the CMD distribution 
jValcarce & Catelanll2008l) : however, it appears extremely un- 
likely that any such disagreements may somehow be due to 
internal variations in Y, since they seem to affect blue and red 
HB stars alike, as also noted by Valcarce & Catelan. 

In Figure[3]we compare our CMD data to ZAHB sequences 
for different Y values in order to obtain a more quantita- 
tive limit on any increase in Y between the red and blue 
HB populations. These ZAHB sequences are all for a fixed 
Z = 0.002 and have been shifted by Amag = -0.16 to match 
the lVandenBergetalJ (l2006h ZAHB for Z = 8.45 x lO"'*, as 
noted previously. 

Clearly, an enhancement in Y by more than 0.01 in the blue 
HB stars would hardly be compatible with the data. In par- 
ticular, we find many more stars below the blue ZAHB for 
F = 0.24 than below the red ZAHB for 7 = 0.23, again with the 
possible exception of the (more uncertain) M,,, {b-y){) plane. 
This strongly suggests that the level of He enhancement is 
most likely less than 0.01 among the cool blue HB stars in 
M3. 

Recall, from Figure 4 in Caloi & D'Antonal (120081) and Ta- 
ble 1 in D'Antona & Caloi (2008), that the bulk of the blue 
HB stars in M3 are predicted, in the helium-enhancement 
scenario, to be enriched in the range between 0.02 to 0.025. 
Such a level of helium enhancement is clearly ruled out by our 
data. Note also that, in the helium enhancement scenario, one 
should expect to see an increase in Y towards bluer colors, but 
this is not supported by our data. 

Finally, Figure |4] compares the empirical and predicted po- 
sitions of blue HB stars in M3 in the \ogg, log Teff plane, for 
three different helium abundances, ranging from 0.23 (panel 
a) to 0.33 (panel c). To produce this plot, we restrict our- 



selves to temperatures cooler than 11,500 K, to avoid the 
well-known complications due to radiative levitation fo r hot- 
ter HB stai-s jGrundahl et al.|[T999HMic"haud et al.ll2008i) . Our 
CMD's indicate that less than 5% of all of the HB stars in 
M3 are hotter than this limit, which corresponds to a color 
(Z7-y)o« -0.025. 

Similarly to what was found in our CMD analysis, the em- 
pirical gravities also seem entirely consistent with a uniform 
value of Y among the blue HB stars in M3. In addition, there is 
no indication of an increase in Y beyond the canonical value, 
a conclusion which becomes even stronger when the data are 
compared with the a-enhanced ZAHB by IVandenBerg et al.l 
(2006) for a y = 0.236, Z = 8.45 x lO'^. Interestingly, the 
available data also suggest that at least the redder blue HB 
stars in Ml 3 (NGC 6205) have a similar helium abundance as 
in M3. The sample size remains relatively small though, and 
therefore data for more stars in both clusters would certainly 
prove of interest to derive more conclusive results. 

5. CONCLUSIONS 

In this Letter, we have shown that high-precision, well- 
populated empirical CMD's, along with spectroscopic grav- 
ities, can be used to pose strong constraints on the presence 
of helium-enriched populations in GC's. Our results strongly 
suggest that any populations that may have formed after a 
main initial burst in M3 likely preserved closely the same he- 
lium content as in the cluster's primordial gas, with a level 
of helium enhancement most likely not higher than 0.01. In 
future papers, we plan to apply similar tests to several other 
GC's. 
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